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Abstract 

We study the dynamics of intracellular calcium oscillations in the presence of proteins that bind 
calcium on multiple sites and that are generally believed to act as passive calcium buffers in cells. 
We find that multisite calcium-binding proteins set a sharp threshold for calcium oscillations. Even 
with high concentrations of calcium-binding proteins, internal noise, which shows up spontaneously 
in cells in the process of calcium wave formation, can lead to self-oscillations. This produces oscilla- 
tory behaviors strikingly similar to those observed in real cells. In addition, for given intracellular 
concentrations of both calcium and calcium-binding proteins the regularity of these oscillations 
changes and reaches a maximum as a function noise variance, and the overall system dynamics 
displays stochastic coherence. We conclude that calcium-binding proteins may have an important 
and active role in cellular communication. 

PACS numbers: 87.17.Aa, 87.18.Tt, 87.19.ln, 05.40.-a 
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Animal cells communicate through the exchange of molecules secreted in the extracellular 
environment and/or by cell-cell contacts through the chemical interactions of molecules 
expressed at the cell surface. Cell communication is at the basis of animal physiology and 
pathology, and the information stored in the environment must ultimately be decoded in the 
cell and propagate to the nucleus where target genes can either switch on or off in response 
to these stimuli. Calcium (Ca^"*") is the most versatile second messenger in living cells 
and translates the information stored in the extracellular environment in time- dependent 
variations of Ca^'*' intracellular concentration. These may take the form of waves, bursts 
and oscillations that propagate in time and space in the cell and through adiacent cells 

Ca^"*" oscillations occur in a large number of cell types such as excitable (e.g. neurons, 
cardiac cells) and non excitable (e.g. hepatocytes, endothelial cells) cells, either sponta- 
neously or after stimulation by hormones, cytokines and neurotransmitters, and they drive 
important functions such as brain and cardiac activity, immune cell activation, hormone 
secretion and cell death jl, 2|. 

Two aspects of intracellular Ca^^ oscillation have received little attention, with some 
notable exceptions : Ca^^ oscillations propagate in the cell in the presence of Ca^^- 
binding proteins that act as buffers, taking up to 99% of Ca^"*" in the cell; the mechanisms of 
Ca^+ wave generation is intrinsically noisy. Here we address both aspects from a dynamical 
perspective. 

Ca^+ concentration within cells is controlled by reversible binding to specific classes of 
proteins that act as Ca^^ sensors and decode the information carried by FM and AM mod- 
ulation of the Ca^^ oscillations [6|]. Many intracellular proteins can bind Ca^^ on multiple 
sites. For example, the activity of key enzymes such as CaM-kinase II is modulated either 
w direct binding to Ca^+ or indirectly by chemical interaction with Ca^^-binding proteins 
21]. Ca^"*" binding to proteins on multiple sites is not driven by enzymatic mechanisms and 
therefore its chemistry obeys to the law of mass action. In addition, the binding sites on 
Ca^+-binding proteins appear to be all equivalent in presence of the EF hand motif which 
allows the specific and reversible binding of Ca^^ ions 

In this paper we investigate Ca^"*" dynamics in cells in the presence of Ca^^-binding proteins 
and noise. To this end, the vitamin D-dependent protein calbindin-D29K (CaL) has been 
chosen as an example of multisite Ca^~^ binding protein, representing a class of several Ca^^ 
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binding proteins that share many biochemical traits with CaL [2j. CaL is a protein with a 
molecular weight of 30 kDa, that is expressed in the cytosol of specific cells of the central 
nervous system such as hippocampal granule cells and cerebellar Purkinje cells 2(]. CaL is a 
member of the EF hand superfamily Ca^'^'-binding proteins. To date, CaL has not yet been 
found to play a role in Ca^+-dependent regulation of enzyme activity although it is expected 
to act as an intracellular Ca^^ buffer. This seemingly minor role may nonetheless have a 
profound effect on cell physiology, and indeed overexpression of CaL in cultured hippocam- 
pal pyramidal neurons affects synaptic plasticity and suppressed post-tetanic potentiation 
10l |. CaL has four reversible binding sites for Ca^"*" and, importantly, the on-off rates for 
the binding reaction have been determined experimentally. Experiments have also shown 
the presence of binding sites with two different affinities and a ratio of sites with high and 



low affinity of 3 to 1 or 2 to 2 

We describe the dynamics of Ca^"*" reversible binding to CaL with the following differential 
system, which is formally the same that we used in a previous study of multisite protein 

nn n 

modification ([7|, 18|; see also [9| for a derivation of these equations): 

^i^^ = -AKnl[CaLo][Ca'^] + Kff^[CaL^] 
d[CaLi] ^ _koffi[CaLi]-3Kni[CaLi][Ca^+] + 



dt 



4koni[CaLo][Ca^+] + 2koffi[CaL2] 



dt 



2koffi[CaL2] - 2koni[CaL2][Ca'+] + 
3koni[CaL,][Ca^+] + SkoffiiCaL^] 
d[CaLs] .3k,^^^lCaL,]-k,^,[CaL,][Ca'+] + 



dt 



2koni[CaL2][Ca^+] + 4koff2[CaL^] 



dt 



-^koff2[CaLi] + kon2[C aL'iWC a 



3^ = -^n^^^ + /([Ca2+],t) (1) 

Ojit (XL 
n=l 

where square brakets denote molar concentrations of each chemical species, CaL^ with i = 
0,1,..., 4 denotes CaL with i Ca^+ bound ions, and /([Ca^^],t) is a function describing 
the time-dependent oscillations of Ca^"*" in the cell. We take the following values for the 
model parameters [lli: fco„i=1.3 ■ 10^M-is-\ koj fi=2.275s-'\ fcon2=7.7 ■ 10^M-is-\ and 
^o//2=39.501 s^^. The system of equations ([1]) has no analytical solutions and thus it must 



be solved with numerical methods. 

We still have to specify the function /([Ca^+],t); there are many different models of 
intracellular calcium oscillations, but this choice is not critical, and here we take the minimal 
and well-known model based on Ca^^-induced Ca^^ release (CICR) as the basic model of 
Ca^+ oscillations which is sketched in figH] (For recent reviews, see also js, 6] and 
references therein). 

The CICR model must be complemented by a stochastic term because the process of 
Ca^+ wave generation is intrinsically noisy Q, [l^. Evidence of intracellular noise comes 
from direct inspection of sampled time-series in different cell types at the mesoscopic scale 
as well as from experimental work which shows that, in the microscopic domain, Ca^"^ waves 
originate from discrete random events, called puffs, occurring at specific sites in the cells 
and composed of a small number of Ca^"*" ions. Several puffs cooperate to the formation 
of a supercritical nucleus that initiatiates a Ca^"'' wave. The probabilistic character of 
nucleation introduces variability into the wave period, with a standard deviation which has 
been experimentally estimated to reach values up to 40% |l4|. Thus, internal noise in Ca^^ 
dynamics cannot be neglected, and we introduce a stochastic term (as in 15|, llTI] that leads 
to the following stochastic differential system for the driving function / = /([Ca^^],t): 

/ = = vo + vi(3 -V2 + V3 + kf[Y] 

^t^] - v,-v,-kf[Y] (2) 



where 



and 



dt 

[Ca2+]" 



f 2 = 



M2 



V3 = V} 



M3" 



and ^(t) is Gaussian white noise with zero mean and variance D. 

In these equations (see also figH]), [Ca^+] is the cytosolic Ca^+ concentration, whereas 
\Y] denotes Ca^+ concentration in the IP3-insensitive intracellular store, vq = 10~^ M s~^ 
is the input rate of Ca^'^ from the extracellular medium, vi = 7.3 ■ 10~^Ms~^ is the pa- 
rameter related to the IP3-modulated release of Ca^"^ from the IP3-sensitive store. Vm2 = 
65 ■ 10"^ M s~^ denotes the maximum rate of Ca^ pumping into the IP3-insensitive store. 
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whereas 14^3 = 500 ■ 10"^ M s"^ is the maximum rate of release of Ca^"*" from that store 
into the cytosol in a process activated by cytosohc Ca^^; K2 = 10~^M, Kji = 2 • 10~^M 
and Ka = 0.9 ■ 10~^M are threshold constants for pumping, release and activation, respec- 
tively; kf = 1 s^^ is a rate constant that regulates the passive, linear leak of Ca^^ from the 
IP3-insensitive store into the cytosol; k = 10 regulates the assumed linear transport of 
cytosolic Ca^"*" into the extracellular medium. The exponents n = m = 2 and p = 4 denote 
the Hill coefficients characterizing these processes. 

The parameter f3 regulates the saturation of the IPS receptor and acts as the control 
parameter which sets the level of the stimulus and varies from to 1. We also introduce a 
parameter r > 1 s that scales the rates Vq, Vi, Vm2, Vms, kf and k. In this way we can tune 
the Ca^^ period of the CICR oscillator to match oscillations observed in real experiments. 

We use equations ([T]) and ([2]) to investigate numerically the dynamic interplay between 
Ca^+, Ca^"'"-binding proteins and noise in the cell. We integrate the stochastic differential 
system with the Euler-Maruyama algorithm [l^, and for stability reasons we take a short 
integration time step At = 0.001s. 

Before studying the influence of internal noise, we investigate the deterministic system's 
dynamics both in the absence and in the presence of CaL. Simulations show that in the 
absence of CaL, the system undergoes two Hopf bifurcations at /5 ~ 0.29 and /5 ~ 0.77. 
Increasing concentrations of CaL damp the Ca^^ oscillations, until they abruptly switch off 
at the critical concentration [CaL]{0)/[Ca'^~^]{0) ^ 1.6 (with f3 = 0.3, see figl2]). Beyond 
this critical CaL concentration no Ca^+ oscillations are observed in the deterministic case. 
However, fluctuations of Ca^+ concentration due to internal noise allow the system to cross 
the threshold again and oscillate even in the presence of supercritical CaL concentrations 
(figj3]). One feature of the simulation outputs in figj3]is that they are strikingly similar to 



time series sampled in real cells 



181 ]. In these experiments, Ca^^ spikes have been observed 



to occur randomly, and the average interspike interval (T) measured independently for hun- 
dreds of cells is correlated to the standard deviation cr in a rather broad range of (T) values 
jisl ]. We find the same result in our simulations (figJSj). 

Thus, both experiments and the numerical results shown in figEj suggest that noise- 
induced Ca^+ spiking and spiking periodicity depend on noise variance, and moreover the 
observed correlation between noise variance and oscillation frequency also suggest that what 

n n 

we are witnessing here is a form of stochastic coherence [17[. As in [17[ we study this aspect 



5 



using the regularity parameter R, which is the ratio between the average and standard 
deviation of the interspike interval (the "period" of the oscillations): 



As already noted in [l7|, the reciprocal of this quantity is just the "coefficient of variation" 
often used in neuroscience as an estimator of the regularity interspike intervals. In our 
simulations (see figS]) we find that regularity peaks at a given noise variance, and that the 
position of the peak also depends on the CaL concentration. Thus the whole system dis- 
plays stochastic coherence, and proteins that bind Ca^+ ions on multiple sites can tune Ca^+ 
spiking in cells and may ultimately regulate cell communication. Remarkably, experiments 
with engineered knock-out mice for CaL expression show severe impairment of motor coor- 
dination suggesting that CaL has an important role in cerebellar functions and intercellular 
communication [19]. We conclude that proteins, such as CaL, may not simply act as passive 
Ca^+ buffers in cells, but be prime actors in the complex play of cellular communication. 
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FIG. 1: Schematic representation of the CICR model. A signal S acts on receptor R and triggers the 
production of IP3 that stimulates the release of Ca^"'' ions from IP3-sensitive intracellular stores. 
IPS regulates the constant flow of Ca^"^ into the cytosol {vif3). Cytosolic Ca^+ions are pumped into 
an IP3-insensitive store (^'2); Ca?'^ in this store (Y) is transported back to the cytosol in a process 
activated by Ca'^"^ itself (fs). Rates vq, k and kf denote replenishment of the IP3-sensitive store 
with extracellular Ca'^"*", the efflux of cytosolic Ca^+ from the cell and passive leak of Ca^+ from 
the store Y into the cytosol. More details can be found in Q]. This classical scheme is modified 
here to take into account the action of Ca^^ protein buffers such as calbindin (CaL) that binds 
Ca^+ ions on multiple sites. Binding follows the law of mass action and is described by rates kon 
and koff. Finally, in our model cytosolic Ca^"*" fluctuates because of internal noise. 



8 



1. 410* , 




1 10: 



+ 

^ ilO:^ 



l5 I I I I I I 

13 as 1 15 2 2.S 3 

FIG. 2: Upper panel: Bifurcation diagram for the deterministic dynamics without CaL: in the 
absence of CaL the system displays two Hopf bifurcations as a function of the control parameter 
p. Lower panel: here /3=0.3, which brings the system beyond the first Hopf bifurcation, into the 
oscillatory regime, however if CaL rises above a critical concentration ratio [CaL]{0)/[Ca]{0) pa 1.6, 
it blocks all oscillations. Thus, CaL sets a sharp threshold for Ca^+ spiking. 
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FIG. 3: Simulation outputs that show Ca^^ spiking in presence of the Ca^~''-binding protein CaL 
and noise. In all panels the initial conditions are: [Ca^+]{0) = 2 • lO"'^ M, [CaL]{0)/[Ca^+]{0) = 2, 



/3 = 0.3. All the other parameters are as described in the text. Top panels: logj^Q D 



-16 M2. 



Middle panels: log^o = -15.4 M^. Bottom panels: log^o = -13.6 M^. Time series are shown 
on the left whereas the plots on the right show the correlation between average interspike interval 
(T) and standard deviation a. These simulation results are strikingly similar to those found in 



actual experiments 
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FIG. 4: Plot of the regularity parameter R vs. noise variance for the indicated [CaL](0)/[Ca^+](0) 
ratios. Since both Ca^^ and CaL concentrations vary greatly among different cell types, we use the 
ratios [CaL]/[Ca'^^] rather than concentrations. The presence of CaL deeply influences the nature 
of the calcium oscillator and leads to a stabilization of interspike intervals: this regularization of 
the Ca^"*" oscillations depends on noise variance and reaches a maximum for a nonvanishing value 
of the noise variance. 
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